ELSEVIER  Journal  of  Power  Sources  88  (2000)  243-249 


www.elsevier.com /locate /jpowsour 


A  comprehensive  model  for  cyclic  voltammetric  study 
of  intercalation /de-intercalation  process  incorporating  charge  transfer, 

ion  transport  and  thin  layer  phenomena 

M.  Noel a’  * ,  V.  Rajendran  b 

a  Central  Electrochemical  Research  Institute,  Karaikudi  630006,  India 
b  Alagappa  Chettiar  College  of  Engineering  and  Technology,  Karaikudi  630004,  India 

Received  29  September  1999;  accepted  17  December  1999 


Abstract 

A  comprehensive  model  involving  intercalation,  competitive  background  reaction  and  de-intercalation  is  described.  The  interaction 
step  is  modelled  as  simple  irreversible  charge  transfer  and  charge  transfer  coupled  with  diffusion.  The  de-intercalation  step  is  modeled  as 
a  thin-film  process.  Overall  cyclic  voltammetric  (CV)  responses  under  all  these  limiting  conditions  are  also  presented  and  discussed.  The 
method  is  also  used  for  comparing  typical  experimental  CV  responses  reported  in  the  literature.  The  method  is  found  to  be  useful  in 
comparing  the  intercalation  processes  under  different  experimental  conditions  and  evaluating  model  parameters.  ©  2000  Elsevier  Science 
S.A.  All  rights  reserved. 
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1.  Introduction 

In  studies  of  the  electrochemical  intercalation  of  ionic 
species  on  graphite  [1]  and  other  electrodes  in  aqueous 
[2-5]  as  well  as  non-aqueous  media  [6,7],  cyclic  voltam¬ 
metry  (CV)  has  been  commonly  used  for  a  quick  estima¬ 
tion  of  intercalation/de-intercalation  efficiency  (IDE)  un¬ 
der  various  operating  conditions.  This  trend  continues  in 
recent  research  work  in  this  laboratory  [8-11],  and  else¬ 
where  [12-14].  Some  reviews  in  this  area  also  highlight 
the  use  of  CV  as  an  important  analytical  tool  along  with 
other  methods  [15-17]. 

Despite  such  extensive  use  of  CV  in  intercalation/de- 
intercalation  studies,  a  comprehensive  theoretical  model 
applicable  for  CV  has  not  been  available  to  date.  This  is 
mainly  due  to  the  involvement  of  different  phenomenologi¬ 
cal  components  in  the  description  of  this  process.  The 
intercalation  process  itself  may  be  charge-transfer  or 
mass-transport  controlled.  Invariably,  the  intercalation  pro- 
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cess  occurs  along  with  a  competitive  parallel  process. 
De-intercalation  on  the  other  hand  is  predominantly  a 
thin-film  surface  process.  Building  up  the  model  for  such 
an  involved  overall  process  would  indeed  be  interesting 
and  challenging.  Recently,  Levi  and  Aurbach  [18,19]  have 
achieved  an  ingenious  simplification  of  the  overall  process 
by  working  at  very  slow  sweep  rates  in  the  range  of  piV 
s_1  [18,19].  Under  these  conditions,  the  intercalation  as 
well  as  the  de-intercalation  process  can  be  modelled  as  an 
overall  thin-film  process,  as  was  developed  earlier  [20,21]. 

The  above  approach  of  working  at  very  slow  sweep 
rates  or  near  quasi-equilibrium  conditions,  however,  has  its 
own  limitations.  Each  measurement  under  these  conditions 
of  course  involves  significantly  longer  time  scales.  Apart 
from  this,  it  is  not  possible  to  look  at  faster  time-dependent 
processes  such  as  charge-transfer  kinetics  and  diffusion. 
For  the  calculation  of  diffusion  coefficients  of  intercalants, 
for  example,  it  becomes  necessary  to  use  other  techniques 
such  as  AC  impedance  [22,23].  Hence,  it  is  worthwhile  to 
develop  a  model  on  conventional  time  scales  involving 
charge-transfer,  counter-ion  transport  as  well  as  thin  film 
kinetics.  This  is  the  objective  of  the  present  work. 
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2.  Model 

Simple  reversible  anionic  or  cationic  intercalation  pro¬ 
cess  on  a  typical  graphite  electrode  may  be  represented  by 
Eqs.  (1)  and  (2),  respectively,  i.e., 

C,  +  C104-^c:cl04-  +  e-  (1) 

C  v  +  Li++  e-^  C“Li+.  (2) 

As  a  general  case  involving  any  lattice,  L  and  any 
monovalent  cationic  species  M+,  the  cationic  intercalation 
process  may  be  represented  as  follows: 

L  +  M++e~^L“M+.  (3) 

The  rate  expressions  are  considered  for  this  overall 
cationic  intercalation  process,  with  the  cationic  intercala¬ 
tion  current  as  positive.  The  extension  of  these  expressions 
for  anionic  intercalation,  however,  is  obvious  and  straight¬ 
forward.  The  intercalation  reaction,  Eq.  (3),  may  be  deter¬ 
mined  by  the  following  processes: 

(i)  a  simple  activated  irreversible  charge-transfer  pro¬ 
cess; 

(ii)  a  reversible  charge-transfer  process  coupled  with  the 
rate  of  diffusion  of  cationic  species  into  the  lattice 
space;  and 

(iii)  an  irreversible  charge-transfer  coupled  with  the  rate 
of  diffusion  of  intercalant  into  the  lattice. 

Apart  from  the  intercalation  process  at  extreme  negative 
potentials,  additional  cathodic  processes  such  as  cation, 
solvent  or  impurity  reduction  and  hydrogen  evolution  may 
occur.  The  overall  IDE  will  decrease  due  to  such  competi¬ 
tive  processes. 

The  intercalated  film  can  be  modelled  essentially  as  a 
thin  film.  The  de-intercalation  process  may,  thus,  be  mod¬ 
elled  as  a  thin  film  process  as  described  by  Levi  and 
Aurbach  [18,19].  In  what  follows,  the  overall  voltammetric 
responses  for  the  intercalation /de-intercalation  process  in¬ 
volving  an  intercalation  current  density  zin,  a  competitive 
background  cathodic  current  zb,  and  a  de-intercalation 
current  zdi  are  involved.  The  intercalation  current  zin  de¬ 
pends  on  the  model  described  above,  i.e.,  processes  (i)  to 
(iii).  The  overall  responses  are  described  below. 

2.1.  Intercalation  rate  controlled  by  charge  transfer  (case 

1) 

In  this  case,  the  intercalation  process  is  defined  by  a 
simple  irreversible  approximation  of  the  Butler- Volmer, 
i.e., 

hn  =  At  =  AiexP[  - «/“«al(  E  ~  ^th)]  •  (4) 

In  this  expression:  ict  is  the  intercalation  current  density 
due  to  charge  transfer;  zol  is  the  exchange  current  density 


for  the  intercalation  process;  Elh  is  the  equilibrium 
(threshold)  potential  for  intercalation;  anal  is  the  corre¬ 
sponding  transfer  coefficient;  /  is  F/RT\  and  E  is  the 
applied  potential. 

The  competitive  background  current  ib  may  be  repre¬ 
sented  by: 

h  =  A2exp  [~nana2f(  E-Eh)],  (5) 

where:  io2  and  ana2  are  the  exchange  current  density  and 
transfer  coefficient  of  the  competitive  cathodic  process, 
respectively,  and  Eh  is  the  equilibrium  potential  of  the 
competitive  process. 

The  relative  contribution  of  the  overall  charge  from 
these  two  processes  ultimately  determines  the  efficiency  of 
overall  intercalation.  The  IDE  may  be  easily  obtained  by 

IDE  =  qin/(qin  +  qh).  (6) 

In  this  expression,  qin  and  qh  values  may  be  obtained 
by  summing  the  z'in  and  z'b  values  at  different  potentials  or 
time  intervals  in  the  voltammetric  method  using  the  trape¬ 
zoidal  rule,  the  charge  expressions  in  accordance  with  the 
above  rule  are: 

4in  =  V2[0’in,l  +  Hn  ,n)  +  2(*m,2  +  hn,3  +  •  •  •  An.(n-l))] 

(7) 

q b  =  V2[(  A,1  +  K,n)  +2(A,2  +  A, 3+  •  •  •  h,(n-  1))]  •  (8) 

In  this  expression,  the  time  interval  ‘h’  for  CV  is  given  by 
Eq.  (9),  where  v  =  sweep  rate  and  n  is  the  total  number  of 
interactions. 

h  =  E/n  =  ut/n  (9) 

The  overall  cathodic  current  zc  at  any  potential  is 
obviously  given  by: 

A  =  hn  +  A  =  At  +  A-  (10) 

The  de-intercalation  current  for  the  thin  film  model  is 
given  by: 

Ai  =  4in  (dd/dt),  (11) 

where  6  is  the  total  mole  fraction  of  intercalation  species 
in  the  film  with  an  initial  value  of  unity.  The  dO/dt  may 
be  described  in  terms  of  a  charge-transfer  rate  constant  ks, 
an  equilibrium  potential  for  de-intercalation  Ed,  (which 
may  or  may  not  be  equal  to  Eth,  see  below),  and  a  lateral 
interaction  parameter  g,  i.e., 

dO/d  t 

=  -ks6exp[-ctn^f((E- Ed)  +  g8/2)] 

+  K0-  0)exp[  an,J((  E  -  £d)  +  gO/2)] ,  (12) 

where  a  na3  is  the  transfer  coefficient  for  the  de-intercala¬ 
tion  process. 
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For  the  reverse  anodic  peak,  the  potential  and  current 
are  defined  by: 

E  =  2  vtm  —  vt  (13) 

*a  =  *b  +  *de-  (14) 

The  effect  of  intercalation  potential  and  the  overall 
voltammetric  response  are  shown  in  Fig.  1.  For  a  fixed 
value  of  background  potential  Eh ,  as  the  difference  (  Elh  — 
Eh)  increases  the  intercalation  charge,  qin,  as  well  as  IDE 
increases  (Fig.  1). 

For  a  constant  difference  in  the  equilibrium  potentials 
(Eth  —  Eh),  the  relative  values  of  exchange  current  density 
of  intercalation,  iol,  and  background  exchange  current 
density,  io2,  may  also  have  significant  influence  on  the 
voltammetric  response.  Obviously,  higher  values  of  iol 
lead  to  higher  intercalation  efficiency. 

The  de-intercalation  potential,  Ed,  is  a  very  important 
process  parameter.  Under  ideal  conditions,  where  intercala¬ 
tion/de-intercalation  proceeds  on  a  host  lattice  without 
any  surface  transformations,  Eth  must  be  equal  to  Ed.  In 
real  situations,  however,  there  is  a  scope  to  consider  the 
possibility  of  different  de-intercalation  potentials.  This  pa¬ 
rameter  influences  the  peak  potential  alone  and  not  peak 
current. 

The  rate  constant  of  the  de-intercalation  process,  ks,  can 
also  influence  the  de-intercalation  peak  in  similar  fashion, 
as  shown  in  Fig.  2.  The  de-intercalation  peak  current  value 
once  again  remains  unaffected  (Fig.  2)  for  slow  surface 
processes. 

The  lateral  interaction  parameter,  g,  can  be  considered 
as  an  overall  lump  parameter  covering  the  interaction 
between  the  intercalated  species,  as  well  as  lattice  transfor¬ 
mations  and  lattice  intercalant  interactions.  This  parameter 
influences  the  de-intercalation  peak  current  and  the  half¬ 
peak  width  of  the  de-intercalation  peak  (Fig.  3). 


Fig.  2.  Theoretical  CV  for  intercalation/de-intercalation  process.  Influ¬ 
ence  of  ks:( a)  ks  =  0.1,  (b)  ks  =  0.01,  (c)  ks  =  0.001. 


2.2.  Intercalation  rate  controlled  by  reversible  charge- 
transfer  and  diffusion  ( case  2) 


Under  the  operating  conditions  of  intercalated  battery 
systems  over  extended  periods,  the  overall  rate  of  interca¬ 
lation  ultimately  depends  on  the  rate  at  which  the  inter- 
calants  diffuse  into  the  host  lattice.  In  most  cases,  even 
under  practical  CV  time  scales,  the  diffusion  can  become 
the  rate-controlling  process.  The  two  stages  of  the  general¬ 
ized  intercalation  reaction  3  may  be  represented  as  elec¬ 
tron-transfer  (step  3a)  followed  by  diffusion  of  cationic 
intercalant  into  the  lattice  (step  3b),  where  L  and  L“ 
represent  neutral  and  charge  lattice  sites. 


U  +  e — U 


L"+M+ 


Diffusion 

=  L"M+ 


(3a) 

(3b) 


The  electron-transfer  coupled  with  diffusion  process  is 
quite  similar  to  the  solution  phase  electron-transfer/diffu- 
sion  problem  treated  earlier  [24],  but  with  one  significant 
difference.  The  concentration  gradient  is  now  on  the  elec¬ 
trode  side  of  the  electrode-electrolyte  interface  rather  than 


Fig.  1.  Theoretical  CV  for  intercalation/de-intercalation  process;  Influ-  Fig.  3.  Theoretical  CV  for  intercalation/de-intercalation  process.  Influ¬ 
ence  of  (is^  —  Eb):(Eth  —  Eb)  =  (a)  100  mV,  (b)  150  mV,  (c)  200  mV.  ence  of  g:  (a)  g  =  —  5,  (b)  g  =  —  3,  (c)  g  =  —  1. 
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on  the  solution  side.  This  problem  is  also  similar  to  the 
electrodeposition  of  metals  like  cadmium  on  a  mercury 
electrode  under  semi-infinite  linear  diffusion  conditions. 
The  concentration  gradient  in  this  process  during  the  an¬ 
odic  dissolution  step  occurs  in  the  interfacial  layer  of  the 
mercury  film  in  contact  with  the  electrolyte  solution.  In  the 
intercalation/de-intercalation  process,  however,  the  con¬ 
centration  gradient  during  intercalation  as  well  as  de-inter- 
calation  steps  occurs  at  the  electrode  surface  region. 

Despite  these  differences,  mathematically,  the  model  for 
all  these  processes,  which  involve  electron  transfer  coupled 
with  semi-infinite  linear  diffusion,  is  the  same.  Hence,  the 
overall  diffusion-limited  current  expression  developed  ear¬ 
lier  [24]  can  also  be  employed  here.  The  overall  intercala¬ 
tion  current,  zin,  for  a  reversible  charge-transfer  coupled 
with  diffusion  (zrd)  is  given  by  [24]: 


hn  =  h,d  =  nFACjDnfvTTx(at).  (15) 

In  this  expression,  C0  refers  to  the  initial  bulk  concentra¬ 
tion  of  the  intercalant  whose  diffusion  towards  the  lattice 
sites  (L“)  controls  the  limiting  rate.  The  above  expression 
was  originally  developed  for  diffusion  of  electroactive 
species  in  a  solution  whose  concentration  is  very  low 
(mM).  With  intercalation  processes,  the  intercalant  concen¬ 
tration  in  the  solution  is  significantly  higher  and,  thus, 
diffusion  of  counter  ions  into  the  host  lattice  from  the 
interphase  is  the  rate-determining  step.  For  this  phenom¬ 
ena,  the  same  rate  expressions  can  be  employed.  For 
example,  this  procedure  is  adopted  in  describing  the  CV 
behaviour  of  metal  ion  deposition  and  subsequent  dissolu¬ 
tion  from  a  hanging  mercury  drop  electrode  [24]. 

In  this  equation,  the  current  function  ^(zzt)  can  be 
expressed  as  an  infinite  series  solution  as  proposed  by 
Nicholson  and  Shain  [24],  i.e., 


AuO) 


oo 


E  (  -  !)i+  V;'exp[  -jnf(  E  -  £th)] . 

7=1 


(16) 


In  a  recent  work  by  one  of  the  authors  [25],  the  proce¬ 
dure  for  summing  up  this  series  by  the  Pade  approximation 
method  using  continued  fraction  has  been  described  in 
detail  [25].  The  Pade  coefficients  for  [10/11]  px,  p2  . . .  pl0 
and  qx,  Pade  q2...qn  are  presented  in  Table  1.  Using 


Fig.  4.  Theoretical  CV  for  intercalation/de-intercalation  process  (reversi¬ 
ble  charge  transfer  coupled  with  diffusion).  Influence  of  D:  (a)  D  =  1 0  ~ 7 , 
(b)  D  =  10“6,  (c)  D  =  5XKT5. 


these  coefficients,  the  function  ^(rz/)  can  be  calculated 
from  the  Pade  expression: 

pxx  +  p2x2  +  ... 
qxx  +  q2x2  +  q3x3  . . . 

where  x  =  (E  —  Eth)n. 

The  calculation  of  the  overall  current  potential  response 
is  once  again  quite  similar  to  the  procedure  described  in 
Section  2.1,  where  the  current  expression  for  the  intercala¬ 
tion  in  Eq.  (4)  is  replaced  by  new  current  expression  for  zin 
(Eq.  15). 

Apart  from  the  parameters  that  influence  z'b  and  zdi,  the 
overall  CV  response  in  this  case  may  be  influenced  by  E{h, 
the  diffusion  coefficient  of  intercalant  in  the  host  lattice 
and  the  sweep  rate.  Typical  CV  responses  showing  the 
effect  of  the  diffusion  coefficient  are  presented  in  Fig.  4. 
The  intercalation  and  de-intercalation  peak  current  and 
peak  potential  values  are  significantly  influenced  by  the 
diffusion-coefficient  and  sweep  rate. 

2.3.  Intercalation  rate  controlled  by  irreversible  charge- 
transfer  and  diffusion  (case  3) 

The  expression  for  the  intercalation  current  [24],  zin,  for 
irreversible  charge-transfer  coupled  with  diffusion  is  given 
by  Eq.  (18).  Based  on  arguments  similar  to  those  proposed 
in  Section  2.2,  this  expression  can  be  employed  for  calcu- 


X(at) 


Table  1 


Pade’s  coefficients  for  calculating  ^(az)  involving  reversible  charge-transfer  coupled  with  diffusion  (case  2) 


i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Pi 

V 

-136.3 

-136.3 

-239.5 

-432.3 

114.9 

-260.3 

138.6 

246.7 

-145.1 

94.8 

56.7 
- 124.3 

217.9 

-255.1 

94.7 

429.4 

11.6 

195.6 

0.25 

30.4 

1.0 
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Table  2 


Pade’s  coefficients  for  calculating  involving  irreversible  charge-transfer  coupled  with  diffusion  (case  3) 


i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Pi 

2829.6 

2274.2 

864.0 

-91.5 

-91.5 

36.4 

55.8 

21.19 

4.1 

0.026 

- 

9i 

1596.4 

4112.7 

4230.7 

1940.0 

121.7 

- 192.0 

30.6 

100.6 

49.5 

11.0 

1.00 

lating  the  counter  ionic  diffusion  coupled  with  irreversible 
charge-transfer  in  to  the  host  lattice. 

hn  =  *r,d  =  nFAC0{a~^DfuTT  x(  bt) ,  (18) 


where:/? T\(bt)  =  £  (  -  1)J+ 1  -  1)! 

7=1 

X  exp [-j'Mnal/((£-  Eth ) 

+  (\/an^)\n(fnrdb  /kh))\.  (19) 

In  this  expression  :b  =  analnfv,  v  =  sweep  rate.  (20) 

The  current  function  x(kt)  can  be  calculated  using  the 
analytical  expression,  which  is  identical  to  Eq.  (17),  for 
which  px  and  q{  values  are  given  in  Table  2  and  v  is 
given  by  Eq.  (21),  i.e., 


v  —  crnal 


E  ~  +  (  RT/an.iF)log{jwD(lb  /£„) 


(21) 


The  calculation  of  overall  current-potential  response  is 
once  again  quite  similar  to  the  procedure  described  in 
Section  2.1,  where  the  expression  for  intercalation,  namely, 
Eq.  (4),  is  replaced  by  the  new  current  expression,  Eq. 
(18). 

Typical  CV  showing  the  effect  of  the  heterogeneous 
rate  constant  on  the  overall  response  are  shown  in  Fig.  5. 
This  parameter  is  found  to  influence  both  intercalation  and 
de-intercalation  peaks  significantly  (Fig.  5).  When  the 
intercalation  process  occurs  quite  close  to  the  background 
process  (for  the  lower  overall  heterogeneous  rate  constant) 


the  intercalation  peak  current  appears  to  be  higher.  This  is 
purely  due,  however,  to  a  merger  of  the  background 
current  along  with  the  intercalation  current. 

Quite  interestingly,  even  slight  changes  in  the  transfer 
coefficient  values  can  influence  significantly  the  overall 
voltammetric  response  (Fig.  6).  The  overall  IDE  can  also 
vary  accordingly.  The  sweep  rate  can  also  influence  the 
voltammetric  response.  At  higher  sweep  rates,  the  interca¬ 
lation  charge  decreases.  By  contrast,  the  overall  IDE  in¬ 
creases  with  sweep  rate. 


3.  Comparison  with  experimental  responses  and  evalu¬ 
ation  of  model  parameters 

A  number  of  systems  show  CV  responses,  which  corre¬ 
spond  to  simple  irreversible  charge-transfer  control.  Typi¬ 
cal  CV  responses  for  fluoride  ion  intercalation  on  a  graphite 
surface  are  shown  in  Fig.  7  [8].  The  model  parameters 
evaluated  from  closely  similar  theoretical  voltammograms 
are  also  presented  for  comparison. 

Experimental  voltammetric  responses  for  fluoride  ion 
intercalation  in  aqueous  (Fig.  7a)  and  aqueous  methanolic 
media  (Fig.  7b)  are  quite  similar.  The  IDE,  however,  is 
shown  experimentally  to  be  higher  in  aqueous  methanolic 
media.  The  voltammetric  simulation  data  (see  caption  of 
Fig.  7)  clearly  indicates  that  the  exchange  current  density 
for  the  background  oxidation  process,  io2,  is  considerably 
lower  in  aqueous  methanolic  media.  The  /01/z'02  value  is 


Fig.  5.  Theoretical  CV  for  intercalation/de-intercalation  process  (irre¬ 
versible  charge-transfer  coupled  with  diffusion).  Influence  of  kh :  (a) 
kh  =  0.001,  (b)  kh  =  0.0001,  (c)  kh  =  0.00001. 


Potential  (mV) 

Fig.  6.  Theoretical  CV  for  intercalation/de-intercalation  process  (irre¬ 
versible  charge-transfer  coupled  with  diffusion).  Influence  of  ana2:  (a) 
« na2  =  0.3,  (b)  ana2  =  0.4,  (c)  ana2  =  0.6. 
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Potential  (mV) 

Fig.  7.  Experimental  (o)  and  theoretical  (-)  CV  for  a  HPC  electrode  in  10 
M  HF:  (a)  aqueous  solution  (/ol  =  0.012  mA,  io2  =  0.01  mA,  Eth  =  1100 
mV,  Eb  =  1200  mV,  ks  =  0.1,  v  =  40  mV,  anal  =  0.4,  g  =  -  2);  (b) 
aqueous  methanolic  (50  vol.%)  solution  (/ol  =  0.005  mA,  io2  =  0.001 
mA,  Eth  =  1100  mV,  Eb  =  1200  mV,  ks  =  0.1,  v  =  40  mV,  anal  =  0.4, 
8  =  ~  2). 

close  to  unity  in  aqueous  media,  but  is  5  in  aqueous 
methanolic  media.  This  lower  background  current  essen¬ 
tially  leads  to  higher  efficiency  in  aqueous  methanolic 
media. 

The  classical  study  of  anion  intercalation  by  Beck  et  al. 
[2]  seems  to  follow  the  model  of  reversible  charge-transfer 
coupled  with  diffusion  rate.  Typical  experimental  results 
are  compared  with  the  models  voltammograms  in  Fig.  8. 
In  strong  acidic  media,  the  anionic  intercalation  process 
proceeds  much  faster  with  considerably  higher  exchange 
current  densities.  The  intercalation  rate  is  significantly 
higher  than  the  diffusion  rate  and  background  oxidation 
rate.  The  threshold  potential  for  CIO/  and  BF4  intercala¬ 
tions  also  appear  to  be  identical.  The  small  differences  in 


Potential  (mV) 


Fig.  8.  Experimental  (o)  and  theoretical  (-)  CV  for  intercalation/de-in- 
tercalation  process  with  CPP  (v  =  1  mV/s):  (a)  4  M  HCIO  (62%) 
(£th  =  1120  mV,  Eb  =  1120  mV,  v  =  1.0  mV,  D  =  0.000005,  Z0  =  0.5, 
ks  =  0.01);  (b)  5  M  HBF  (50%)  (£th  =  1120  mV,  Eb  =  1120  mV,  v  =  1.0 
mV,  D  =  0.000005,  Z0  =  0.5,  ks  =  0.01),  a/7al  =  ana 2  =  ana3  =  0.4. 
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Fig.  9.  Experimental  (o)  and  theoretical  (-)  CV  of  HOPG  in  0.05  M  Mg 
C12-DMSO  solution  (£th  =  -  1200  mV,  Eb=  -  1500  mV,  Zo2  =  0.00005 
mA,  ks  =  0.01,  kb  =  0.000001,  D0  =  0.00005,  a«al  =  ana2  =  ana3  = 
0.4).  (a)  v  =  5  mV /s,  (b)  v  =  8.33  mV /s,  (c)  v  =  16.66  mV /s. 


the  voltammetric  responses  are  due  to  slight  changes  in  the 
diffusion  coefficients  of  the  intercalants  in  the  graphite 
lattice. 

Cationic  intercalation  studies  by  Maeda  and  Touzain 
[12]  can  be  cited  as  an  example  involving  irreversible 
charge-transfer  coupled  with  diffusion  (Fig.  9).  In  this 
case,  both  the  cathodic  and  anodic  peak  potentials  shift 
considerably  with  sweep  rate.  This  shows  that  the  charge- 
transfer  rate  of  intercalation  (z‘ol)  and  de-intercalation(/cs) 
control  the  overall  process  in  addition  to  the  diffusion 
coefficient.  It  is  also  expected  that  the  intercalation  of  a 
bivalent  cation  like  Mg2+  with  its  solvation  shell  would 
indeed  be  quite  difficult  and,  hence,  rate  controlling. 


4.  Conclusions 

A  simple  model  involving  intercalation,  the  competitive 
background  process  and  the  de-intercalation  process  is 
presented  in  this  work.  The  computation  procedure  for  the 
overall  voltammetric  response  for  this  process  has  also 
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been  developed.  This  has  become  possible  specifically 
because  of  the  availability  of  simple  analytical  expressions 
for  the  current  functions,  which  involve  diffusion  pro¬ 
cesses  [25].  The  model  is  also  found  to  be  useful  in 
evaluating  the  process  parameters  by  comparing  these  with 
the  experimental  voltammetric  response.  This  procedure 
can  be  easily  employed  for  studying  a  wide  range  of 
intercalation  processes  using  one  of  the  simplest  and  rou¬ 
tine  electroanalytical  techniques,  CV.  Further  extension  of 
this  model  to  cover  multiple  stages  of  intercalation /de-in- 
tercalation  is  also  possible.  There  is,  of  course,  consider¬ 
able  scope  for  further  work  in  this  direction. 
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